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^ ' Abstract 

^ ■ Supergravity backgrounds dual to a class of exactly marginal deformations of A/" = 4 su- 
I persymmetric Yang-Mills can be constructed through an SL{2, M) sequence of T-dualities 
Qh! and coordinate shifts. We apply this transformation to multicenter solutions and derive 
r-| , supergravity backgrounds describing the Coulomb branch of A/ = 1 theories at strong 
>• . 't Hooft coupling as marginal deformations of A/" = 4 Yang-Mills. For concreteness we 
^ ■ concentrate to cases with an SO (A) x 5*0(2) symmetry preserved by continuous distribu- 
' tions of D3-branes on a disc and on a three-dimensional spherical shell. We compute the 
expectation value of the Wilson loop operator and confirm the Coulombic behaviour of 
the heavy quark-antiquark potential in the conformal case. When the vev is turned on we 
find situations where a complete screening of the potential arises, as well as a confining 
regime where a linear or a logarithmic potential prevails depending on the ratio of the 
quark-antiquark separation to the typical vev scale. The spectra of massless excitations 
on these backgrounds are analyzed by turning the associated differential equations into 
Schrodinger problems. We find explicit solutions taking into account the entire tower of 
states related to the reduction of type-IIB supergravity to five dimensions, and hence we 
go beyond the s-wave approximation that has been considered before for the undeformed 
case. Arbitrary values of the deformation parameter give rise to the Heun differential 
equation and the related Inozemtsev integrable system, via a non-standard trigonometric 
limit as we explicitly demonstrate. 
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1 Introduction 



The AdS/CFT correspondence is the first exphcit reahzation of the long suspected de- 
scription of the perturbative expansion of an SU{N) gauge theory in the large limit 
in terms of a certain string theory [1] . As it relates the weak coupling regime of A/" = 4 
supersymmetric Yang-Mills to the strong coupling regime of type IIB string theory com- 
pactified on AdS^ x S^, and vice versa [2, 3, 4], a complete formulation of the corre- 
spondence would require a precise knowledge of the strong coupling limit of each theory. 
However, in spite a detailed proof still remains a challenge, diverse tests of the correspon- 
dence, even beyond the supergravity limit from the study of sectors with large quantum 
numbers, have lead to fascinating insights on both the gauge and gravity sides. 

A natural step forward towards a more complete understanding of the duality is the 
extension of the correspondence to less symmetric theories. An appealing candidate is 
the case of the exactly marginal deformations of A/" = 4 Yang-Mills preserving J\f — 1 
supersymmetry [5], because they contain a continuous deformation parameter. A gravity 
dual was recently derived for the Leigh-Strassler or /3-deformation of the Af — A theory 
through an 5'L(2,R) sequence of T-duality transformations and coordinate shifts [6]. 
The construction of this deformed background has uncovered an interesting set of new 
tests of the correspondence. A three-parameter deformation of AdS^ x S'^ dual to a 
non-supersymmetric marginal deformation of A/" = 4 Yang-Mills was soon after found, 
as well as a Lax pair for the bosonic piece of the string in the case of real deformations 
of the background, thus implying integrability of the corresponding bosonic string sigma 
model [7]. In addition, the energies of semiclassical strings rotating with large angular 
momenta in the deformed background were then compared to anomalous dimensions of 
large gauge theory operators in the /3-deformed — A theory [8] . Diverse features of the 
deformed theories and some related spin offs have been further explored recently, both 
on the gravity [9]- [13] and on the field theory [14]- [20] sides of the correspondence. 

In this paper we will extend the construction in [6] to supergravity backgrounds de- 
scribing the Coulomb branch of marginally deformed J\f — 4 Yang-Mills. We enter the 
Coulomb branch of the — A theory when the SO (6) scalar fields acquire Higgs expec- 
tation values. On the gravity side the non- vanishing scalar expectation values correspond 
to a multicenter distributions of branes. We will use some of these backgrounds to find 
a gravity dual for the Coulomb branch of the /3-deformation oi — A Yang- Mills. The 
plan of the paper is the following: In section 2 we will apply a sequence of T-dualities and 
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coordinate shifts to general supergravity backgrounds with at least a [/(l)x[/(l) x[/(l) 
global symmetry group. In section 3 we consider two solutions corresponding to two 
different brane distributions, with SO{2)xSO{2)xSO{2) and SO{4)xSO{2) global sym- 
metry. Section 4 is devoted to the issue of supersymmetry for the marginally deformed 
backgrounds. In section 5 we probe the geometry of the deformation by computing 
the expectation value of the Wilson loop operator along the transversal space for a dis- 
tribution of D3-branes on a disc and on a three-dimensional spherical shell preserving 
SO {A) xSO{2) symmetry. Section 6 contains a detailed analysis of the spectra of mass- 
less excitations by studying the Laplace equation in the backgrounds with SO{4:)xSO{2) 
symmetry as well as the conformal case. The regime of arbitrary values of the deforma- 
tion parameter leads to the Heun differential equation, which we relate to a generalized 
trigonometric limit of the Inozemtsev integrable system. We conclude in section 7 with 
a discussion and further directions of research. 



2 The general set up 

In this section we will construct supergravity backgrounds dual to the Coulomb branch 
of marginally deformed A/" = 4 supersymmetric Yang-Mills following the SL{2,M.) se- 
quence of T-dualities and coordinate shifts introduced in [6]. We will therefore start 
from supergravity solutions modeling the Coulomb branch of the J\f — 4 theory at strong 
coupling. In these backgrounds only a metric and the self-dual 5-form are turned on. 
The corresponding expressions can be written in terms of a harmonic function in H^, i.e. 
the transverse space to the branes. The metric has the form 

dslo = H-^/^ri^^dx''dx'' + H^/^dxidxi , // = 0, 1, 2, 3 , i = 1, 2, . . . , 6 , (2.1) 

the self- dual 5-form is given by 

Fs = dA4 + *iodA4 , ((i^4)oi23i = -diR-' , (2.2) 

and the dilaton is a constant $o- The harmonic function H is in general given by 

where the density of the brane distribution is normalized to unity and should be positive 
definite. We have already taken the field theory limit so that asymptotically the space 
is AdS^ X with radius — AirgyM^ string units. The distribution of the brane 
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centers breaks the 5*0(6) global symmetry of the background. Generically this breaking 
is complete, but, in the cases we will be interested, a smaller subgroup of the 7^-symmetry 
is retained. In particular in order to perform the combination of dualities in this paper 
and the same time preserve some supersymmetry, we will need at least a group isomorphic 
to U{lf. 

Let us start with a general background in which the ten spacetime coordinates are 
split into a seven-dimensional part parametrized by , with / = 1,2, ... ,7, whereas the 
remaining three coordinates form a 3-torus parametrized by the angles (pi, i = 1,2,3. We 
therefore take the following metric 

3 

dslo = Gij{x)dx^dx-^ + J2^ii^)d(l>t , ^=1,2,..., 7, (2.4) 

1=1 

where, as indicated, the three positive definite functions Zi could depend on the transverse 
to the torus coordinates, the x^'s. In the case that the 0j's parametrize the Cartan torus 

of an undeformed 5-sphere, the z*'s sum up to unity. In all other cases in which the 
five-sphere is deformed they do not obey any restriction except of course those arising 
from preserving supersymmetry and satisfying the field equations. The ranges of the 
angle variables, the associated Killing vectors and their norms are 

0, e(0,27r), 6 = (0,1,0,0), 6 = (0,0, 1,0), C3 = (0,0,0,l) , C-=^. ,(2.5) 

where i — 1,2,3. In general there are degeneration surfaces where at least one of the 
norms vanishes. In our case they are defined by the equations Zi — 0. With the above 
choice for the ranges of the angular coordinates, regularity of the metric and absence 
of conical singularities at the degeneration surfaces requires that the associated "surface 
gravity" (important as a notion in black hole solutions with Lorentzian signature) equals 
to one. In our case this means that 



r\2 : 



zi=o Azi 



Zi=0 



~ 1 , (no sum over i — 1,2,3) . (2.6) 



Note that for simplicity we have not allowed mixing terms. Nevertheless such terms 
can also be included to allow for more general backgrounds. Under the above assumptions 
the self-dual 5-form is given by 

Fs = dC^^^ A d(t)i A #2 A d(t)2 + ^ *7 dC^^^ , (2.7) 

^JZ\Z2Zz 

for some 1-form C*^^-* = C^pdx'^ . This is clearly self-dual and the 1-form C^^-* should 
be such that the exterior derivative of the second term is zero. Hence, there must be a 
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4- form C'-^-' such that 

dC^^) = ^ *7 dC(i) . (2.8) 
Consider now the change of variables 

with inverse 

^1 = ^(01 + 02 - 203) , <^2 = ^(02 + 03 - 20l) , ^3 = ^(01 + 02 + 03) • (2.10) 

This transforms the three-dimensional part of the metric to 
3 

^iO?0i = gijdcpidipj = {z2 + zz)d^p\ + {zi + Z2)dipl + (2:1 + 2^2 + Z3)dipl 

i=l 

+ 2z2d(pidip2 + 2{z2 — Zs)d(pidip3 + 2{z2 — Zi)d(p2d(p3 ■ (2-11) 

The U{1) X U{1) global symmetry is generated by shifts of the angles ipi and ip2- This 
shifts correspond to transformations of the three complex superfields of the N — ^ SYM 
theory as 

($i,$2,$3) ^ (*i,e*"i$2,e-*"i*3) and (e"*"^*!, e*"^$2, $3) , (2.12) 
under which the superpotential of the marginally deformed theory 

W = Tr(e'''^$i$2^3 - e-^''^$i$3$2) , (2.13) 

remains invariant. The 7?.-symmetry will then correspond to shifts of the remaining angle 
</73. Since all fields transform in the same way under this symmetry the superpotential is 
not invariant, and transforms with a weight as it should be. The supercharges preserved 
by the original background fit into representations of the total U{1) x U{1) x U{l)'ji 
symmetry. The solutions that will be generated below by the combination of T-dualities 
and a coordinate shift will preserve a fraction of the original supersymmetry because 
they will not involve the angle ip3 and therefore the U{l)'ji symmetry will not mix with 
the global U{1) x U{1). This will be examined in detail in section 4. 

The vacuum structure of the theory (the Coulomb branch) is described by the conditions 
$^$2 = q^2^i , q = e"^'""^ , and cyclic , (2-14) 
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which are vahd for large (exact for U (N)). For general values of 7 this requires traceless 
N X N matrices, where in each entry at most one of them is non-zero.^ We note that in 
the undeformed theory any tracelessA^ x N matrix is allowed. 



2.1 The deformation 

We will now deform our background with an SL{2, R) C SL{3, M) transformation of the 
complete SL{3, R) x SL{2, R) duality group of type-IIB supergravity compactified on the 
global C/(l) X C/(l) torus. ^ We repeat here essentially the steps of [6] and [7], but we will 
no longer consider the conformal constraint Y2i — i- In order to clarify the derivation 
we will include all details. 

As a first step we will perform a T-duality along the Lpi direction. In general, under T- 
duality the fields in the NS-NS sector form a closed set and transform among themselves. 
Hence for these we may use the standard rules [24]. In contrast, the transformation rules 
for the R-R sector fields involve those in the NS-NS sector (see, for instance, [25]). We 
find that the non-zero components of the metric and antisymmetric tensors are 

1 ~ Z1Z2 + Z1Z3 + Z2Z3 
9ii = — : , 922 = ■ , 

Z2 + ^3 Z2 + ^3 

933 = ^1 H : , 5-23 = : Zi , (2.15) 

Z2 + Z3 Z2 + Z3 

I Z2 7 Z2 — Z3 

bi2 = ■ , bi3 = ■ ■ 

Z2 + Z3 Z2 + Z3 

In addition we obtain for the dilaton 

e-'* = -^2 + -^3 , (2.16) 

and for the R-R 3-form 



^fLs = 3^1'^ ■ (2-17) 



^Studies on the Coulomb branch of the theory from the gauge theory side have been done in [21, 22], 
where, in addition to the generic behaviour we study here with gravity duals, it was shown that there 
exist exceptional Coulomb branches for sp(X'ial values of the deformation parameter. 

^Prototype examples where such transformations can be performed are provided by four-dimensional 
NS-NS backgrounds with two commuting U{1) isometries. Starting with the exact string background 
corresponding to the SU{2)/U{1) x SL{2,JEi)/U{l) coset WZW models we recognize, after the trans- 
formation, the background describing NS5-branes uniformly distributed on a circle in their transverse 
space [23]. A similar case to this was also considered in [6] where one starts with IR^ x IR^ in polar 
coordinates, instead of the coset models. 
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We will now perform a coordinate shift 

^2+ 1V\ ■ (2.18) 

Then we find for the metric (we use upper case letters to denote the tensors after the 
shift) 

^ 2~ n-l 1 n ~ -i.:2 + "'i "':-! + -22:3 

^11 =511+7 5*22 = ■ , G22 = 5*22 = ■ , 

^2 + -2^3 Z2 + Z3 

G33 = 5-33 = ^1 H \ , <^23 = 92Z = \ ^1 , (2.19) 

Z2 + -^3 -22 + Zz 

~ _ / 22:2^3 \ ^ _ ^1^2 + ^1^3 + ^2-^3 

<^13 = 7fl'23 = 7 \ 2^1 , <-^12 = 75'22 = 7- 



Z2^Z'i J Z2 + Z^ 

where for notational convenience we have defined 

= 1 + 7^(zi^2 + -^1-^3 + -22-^3) . (2.20) 

The antisymmetric tensor remains unchanged under the coordinate shift, i.e. Bah = ^ab- 
For the R-R 3-form this shift produces the non-zero components 

<1.3 = ^Cf^ ' 4l\^3 = 37^^?' ■ (2.21) 

Next we perform again a T-duality transformation along the </7i direction and find for 
the metric 

G,,- = G (5,,- + 97'ziZ2^3Mj,3) , (2.22) 
and for the antisymmetric tensor 

-B12 = 7<^(^i^2 + 2:12:3 + 2:22:3) , 

-B13 = 7G'(22;22;3 - Z1Z2 - Z1Z3) , (2.23) 

-B23 = 7<^(2:i^3 + 2:22:3 - 22;i2;2) . 

These expressions for the metric and antisymmetric tensor are written in the (pi coor- 
dinate system. Returning now, with the help of (2.10), to the original coordinates we 
obtain the final result 

" 3 

Zid(f)'f + ^'^ ZiZ2Z^{d(f)i #2 + #3)' 



ds\Q = Gijdx^dx'' G 

.1=1 



(2.24) 



For the antisymmetric tensor we get 

B = '-fG{ziZ2d(j)i A d(j)2 + cyclic in 1, 2, 3) (2.25) 



7 



and for the dilaton 



g2$ ^ g2*o G . (2.26) 
Finally the non- vanishing components of the R-R fields are given by 

A(') = -7C(i) A {d(j)i + d(t>2 + d(t>^) , 

^(4) ^ Q(jii) A d4)i A d02 A #3 + C(^) . (2.27) 

The 5-form is then 

Fg = dA^^^ -dBA ^(2) ^ G'ciC7(i) A d(Pi A ci02 A dc^^ + ^ *7 dC^^^ , (2.28) 

where we have used the identity dG = —j'^[d{GziZ2) + cyclic in 1,2,3]. Its self-duality 
with respect to the deformed metric (2.24) can be readily verified. 

As usual, the supergravity description is valid if the curvature of the metric remains 
small compared to the string scale. Let be the overall scale, in string units, of the 
metric (2.4), which implies that both the Zi^s and the scvcn-dimensional metric scale like 
and that the combination 7/?^ is kept constant. Hence, we require that R ^ 1, but 
in addition wc should make sure that the metric does not degenerate at arbitrary points 
due to the combination 7/?^ becoming large. Indeed, if 'yR^ 3> 1 then we see that the 
term Zjci^f in the metric (2.24) scales like l/(7i?)^ and therefore we should require that 
this scale is large compared to unity.^ To summarize, for the supergravity description to 
be valid at a generic point of the manifold we should have 

> 1 and 7i?2 = 7 < . (2.29) 

This is the same condition as that obtained for the conformal case in [6] where it was 
also noted that the last condition is sufficient for the 2-torus parametrized by ipi^2 to 
stay much larger than the string scale after the T-duahties. This can be shown to be the 
case here as well by noting that after the combined T-duahties and the coordinate shift 
transformation, the volume of the 2-torus is 

R^ 

Vol(2-torus) = G{ziZ2 + Z1Z3 + -^2-^3)^^^ - z , (2.30) 

l + Y 

^In the opposite case it is obvious that curvature invariants, such as the square of the Riemann tensor 

for the deformed metric (2.24), will blow up at a generic point of the manifold and the supergravity 

description will no longer be valid. Note also that string loop corrections are negligible at a generic point 

if e*" <C 1, independent of the deformation. 
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where the last expressions indicates, schematically, the way the volume depends on the 
parameters R and 7 at a generic point. In order for it to be larger than unity, the constant 
7 should be smaller than R. 

The second of the conditions in (2.29) means that 7 <^ for the supergravity approxi- 
mation to be valid. Then the deformation parameter q in (2.14) becomes unity. However, 
this does not mean that we should admit arbitrary vev distributions to parametrize the 
Coulomb branch of the theory, but only those corresponding to scalar fields satisfying 
conditions (2.14). 

We finally note that the periodicities of the angular variables 0^ remain intact in the 
deformed background. To see that, note that the norms of the Killing vectors are now 
given by 

ei^G{z, + ^''z,Z2Z^) , i = 1,2,3. (2.31) 

A simple computation, using the fact that = (for at least one Zj), shows that 
dj^f = djZi H — • , where the dots indicate terms that vanish at the degeneration surfaces 
faster than the indicated first term. Using then the definition in (2.6), we find that indeed 
the "surface gravity" equals one, — 1. 

3 Brane distributions 

In this section we will apply the above construction to derive the explicit form of the 
supergravity background describing the Coulomb branch of marginally deformed H — 

4 supersymmetric Yang-Mills for two different brane distributions, namely a uniform 
continuous distribution of D3-branes on a disc, and on a three-dimensional spherical 
shell. These were first constructed as the extremal hmits of rotating D3-brane solutions 
in [26, 23]. They were also used in several investigations in the literature within the 
AdS/CFT correspondence starting with the works of [27, 28] and belong to the rich 
class of examples representing continuous distributions of M- and string theory branes 
on higher dimensional ellipsoids [29, 30, 31]. 
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3.1 Solutions with 50(2) x 50(2) x 50(2) symmetry 



In this case the transverse space coordinates are parametrized as 



xi\ — 7- . n /cos0i 
■s/r^ — bi sm U 



X2 1 \ sm (pi 



^1 = -x/r^ — 62 cos^sin-^ [ ^ 1 , (3.1) 
J ysin 02 

V?^ — O3 COS ^ COS -0 

y ysm 
where 6^, i = 1, 2, 3 are some real constants. The ranges of variables are 

r > max(6i, 62, ^3) , < V' < | , < 0i,2,3 < 27r . (3.2) 

In this coordinate system the metric is given by 



+ r^H^/^ I Aid9^ + A2 cos^ 9d^'^ + 2 "^ cos ^ sin Q cos sin V'ci^c^V' (3.3) 



.^2 - ^'3 



+ (1 - ^) sin^ dd^\ + (1 - ^) cos^ Q sin^ i,d^\ + (l - ^) cos^ Q cos^ i>d^l^ , 



where the diverse functions are defined as 



H 



A = 1 — cos^ ^ — -I (sin^ 9 sin^ + cos^ ~ "I (sin^ 6 cos^ V + sin^ 

+ ^ sin^ e + ^ cos^ 6 sin^ V' + ^ cos^ 6 cos^ V' , (3.4) 



Ai = 1 ^ cos^ ^ 1 sin^ ^ sin^ -0 1 sin^ ^ cos^ ^ ? 



1 COS ip 5- sm ip 



This metric can be interpreted as the supersymmetric limit of the most general non- 
extremal rotating D3-brane solution [32, 33], with bi, 62 and 63 related to the three 
rotation parameters after an adequate Euclidean continuation. It was also derived as 
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a domain wall solution within five-dimensional gauged supergravity and then uplifted 
to string theory [29]. We have also taken the field theory limit so that the space is 
asymptotically AdS^ x S^. 

Comparing now (3.3) with (2.4) we see that 

^1 = AV^ (^1 - ;:2 j sin ^ ' 

^2 = XV^(^l-;ljcos'^ sin^V', (3.5) 

^3 = ^ (^1 - ^ J cos e cos V ■ 

The marginally deformed SO{2) x SO{2) x SO{2) solution is then obtained by reinstaUing 
these expressions in (2. 24)- (2.28).^ 

Consider now the shift 991,2 — > ^1,2 + cti,2- Then we easily check that the complex 
coordinates 

wi = xi+ 1x2 , W2 = + ix4 , W3 = X5 + ixQ , (3.6) 

transform as the corresponding complex superfields in (2.12), and therefore there is in 
that respect agreement with the deformed gauge theory. 

3.2 Solutions with 5'0(4) x 50(2) symmetry 

These solutions can be obtained by letting bi = —Vq and 62 = &3 = into the various 
general expressions of the previous subsection. The radial variable obeys r > and the 
metric (3.3) now becomes 

ds' = H-'/^^dx^dx-^ + Hi/2^1±lt^ dr' 

-\-rQ 

+ //V2 ((^2 ^ ^2 ^Qg2 q^^^q2 ^ (^2 ^ ^2) g.^2 ^^^2 ^ ^2 ^^^2 q^^2^ ^ (3 7) 

where the harmonic function reduces in this case to 

^ " r2(r2 + r2cos2e) ' ^^'^^ 
and the 3-sphere line element is 

dnl = di)'^ + sin^ dcg + cos^ V #3 ■ (3-9) 



*See also [34] for related work using (3.3). 
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The forms necessary to compute the NS-NS and R-R field strengths are 

+ Tq cos"' 

C7(4) ^ -R-^r^{r^ + rl cos^ 9) dt A dxi A dx2 A cixg . (3.10) 

In our parametrization the configuration corresponds to a set of D3-branes uniformly 
distributed along the X1 — X2 plane (equivalently for r = and ^ = 7r/2) on a disc of 
radius tq. In this case we find 

^ R^(^l + !i^Sfl^ '^'cos^^ sin^V^, (3.11) 

= 1^1 + i^!^^ '^'cos^e cos^V' . 

The marginally deformed 5*0(4) x 5*0(2) background is obtained once we enter these 
expressions in (2. 24)- (2.28), with the function G given now by 



0-^ = 1 + 7^ , " , ^ \(r^ + rl) sin^ 9 + cos^ 9 sin^ ^ cos'^ ip] . (3.12) 
+ Tq cos^ fc/ '- 

Finally, let us note that for Tq —Tq we get a distribution of branes on the surface 
of a four-sphere in the x^, . . . ,Xq space of radius Tq (equivalently for r — Tq and ^ = 0). 
In that case the radial variable obeys r > Tq. 

As a general comment we note that the background metric is singular in the location of 
the distribution. The reason for this is that in this place the continuum approximation of 
the distribution breaks down and it should be replaced by its discrete version. If we place 
at each center N^cntev D3-branes such that 1 ^ iVccntcr ^ N then the gravity description 
of the gauge theory is still valid with the background corresponding to AdS^ x in the 
J\f — 4 case, and its exactly marginal deformation in the J\f = 1 case. The radius of the 
space in this case is much smaller that in the continuum case. It is given by 47rg'YM-^center 
in string units, but nevertheless it can still be taken to be much larger than the string 
scale. 

In the remaining part of this article wc will restrict to the simpler case of solutions with 
5*0(4) X 5*0(2) symmetry, or even to the conformal case. The latter follows immediately 
when the vev parameter tq = 0. Then the various expressions for the background fields 
reduce to those in [6]. 



■'0 
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4 Supersymmetry 



In this section we investigate the issue of supersymmetry for the deformed supergravity 
backgrounds and explain in detail the origin of their reduced supersymmetry. We start 
by explicitly showing that the solution for the Killing spinor in the undeformed case can 
split into a part which is a singlet of the U{1) rotations corresponding to the angles (fi 
and (p2: and a part orthogonal to that. After the T-dualities and the shift only this part 
survives and remains a Killing spinor of the deformed theory. For any multicenter metric 
of the form (2.1) the Killing spinor is of the form (this is a mere consequence of the 
supersymmetry algebra. See, for instance, [35]) 

e = H-'/\ , (4.1) 

where eo is a constant spinor subject to the projection 

= eo , (4.2) 

where the indexes refer to the directions along the brane. It is important to realize that 
for a different form of the metric such as, in the case of our examples, (3.3) and (3.7), 
the constant spinor will acquire a coordinate dependence. This can be found by solving 
the covariant version of the condition diCo — 0, that is 

diSo + ufTabeo^O , i = l,2,...,6, (4.3) 

valid in all coordinate systems. We will concentrate on the background with 5*0(4) x 
SO (2) global symmetry in the disc case. Then we have the flat metric in given by 



ds^e = ^'+^^gcOs'^ ^^2 ^ (^2 ^ ^2 ^^^2 q^^q2 ^ ^^2 ^ ^2^ ^-^2 q^^2 

r + Tq 



+ cos^ e{dip'^ + sin^ ip d(f)l + cos^ ip d(f)l) , (4.4) 
as this is read off from (3.7). Using as a frame basis 



= /?;^jHj_cc^_^ ^ ^2 ^ / ^2 _^ J.2 pQg2 Q d9 , = r cos 61 di/j , 
V + Tq V 

= \Jr'^ -\-rQ sin 6 d(pi , = rcosOsinip d(p2 , = rcos^cos-^ d4>3 , (4.5) 

we find that the spinor eo is given by 
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where constant spinor and where we have defined 

tjcin \ 

(^2 ^^2-) 1/2 J '^i = ^24 , (72 = , (73 = Tie . (4.7) 

We are interested in extracting the part of the spinor that is invariant under variations 
of the angles ipi and </72- After rewriting the spinor eo in the coordinate system and a 
simple computation we find that the required spinor invariant under variations of </?i,2 is 
given by 

-r„. — p5/M)ri2pf i[((Tl+(T2+(T3)¥'3 + (o-2-0-l)¥'2+(o-2-0-3)¥'l]^ 

"^0,inv — 00 cQjinv ) 

^ e^^(^''^)r-etri3ei-3v.3g^.^^ . (4.8) 
The constant spinor eo,inv in terms of eo is given by 

eo,inv = - ^1^2 - (7i(7z - (72(Jz)eQ , (4.9) 
where the pref actor acts as a projector and by construction we have 

'-'^l^O.inv '^2^^0,inv '^S^O.inv ■ 

(4.10) 

The part of the spinor that is not invariant under the variations of ipi^2 will not survive 
the T-dualities [36, 37, 38] because these arc performed with respect to ^pi and as the 
latter has been shifted by a (fi component. This conclusion holds within the supcrgravity 
approximation, but in a string theory context we expect to have the broken part of 
supersymmetry realized with operators having no field theory analog [37]. Related to 
(4.2) and (4.10) projections are expected to arise from a careful examination of the 
Killing spinor equations for the deformed background in the generic case when the vev's 
are turned on. Being three independent conditions on the spinor we are left with Af — 1 
supersymmetry. In the conformal limit (4.2) is not necessary, leading to the J\f — 1 
superconformal case. 



5 Wilson loops and the gg-potential 

In this section we will evaluate the Wilson loop operator in the Coulomb branch of 
marginally deformed Af — A supersymmetric Yang-Mills using the 5*0(4) x SO (2) back- 
ground we have described in section 3.2. According to the prescription in [39, 40], the 
expectation value of a Wilson loop in the field theory can be computed by minimizing 
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the Nambu-Goto action for a fundamental string in a given supergravity background. 
The Wilson loop is constructed by pulling one brane apart from the brane distribution, 
thus giving an expectation value $ to a Higgs field. The quarks are the infinitely massive 
W-bosons connecting the brane distribution to the far away brane. In order to probe the 
geometry on the deformed background, we will introduce a relative angle between the 
quarks by giving expectation values $i and $2 to two C/(l) factors in the global gauge 

— * — * — * 

group. In this way we introduce two relative angles 9i — and the ends of the 

Wilson loop, corresponding to the position of the massive quarks, are located at r = 00 

— * — * 

and two different points 9i and 62 on the transversal space. We will, in particular, choose 
coordinates such that the path joining 9i and 62 is parametrized by 4>2 and 03 in (3.9). 
We will take a trajectory with 

TT 

^ = 0, ip — — , 02 = 03 = 0, 2:2,3 = constant . (5.1) 

This is consistent with the corresponding equations of motion provided that the conserved 
angular momenta for 02 and 03 coincide.^ Then, setting these values in (3.7) the reduced 
four-dimensional metric becomes (in the Euclidean time) 

with 

H = ^2(^2 + ^2) ' C = 1 + 4(^2 + ^2) ■ (5-3) 

We must note that there is also a contribution to this background from the antisymmetric 
tensor field. However, as we will only consider static configurations, it will not contribute 
to the string action. 



5.1 Wilson loops along transversal space 

Before constructing the Wilson loop for the deformed SO {4) x SO {2) solution, we will 

present a discussion which can be used in similar situations when other more general 

supergravity backgrounds are considered. In order to confirm with established notation 

in the literature, and also emphasize that it has the meaning of energy on the field theory 

side, we will use u to denote the radial variable r in all the computations in this section. 

^We may consider a trajectory with 9 = 7r/2 which is also consistent with the equations of motion. 
However this choice is not sensitive to the deformation, since the dependence on the deformation pa- 
rameter 7 drops out. The results for this trajectory in the particular case of no angular separation for 
the quark-antiquark system (/ = below) can be found in [28]. 
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The Nambu-Goto action for a static fundamental string stretching along a great circle 
in the transversal background is 

^ = ^ y dx^/g{u){d.,uf + + h{u){d,(Pf , (5.4) 

where g{u) — grrguu, f{u) — R'^QttQxx and h{u) — gTrg^xj)- The factor T comes from the 
time integration for a static configuration, because we have taken a rectangular Wilson 
loop on the boundary, with one side of length L along the space direction and one of length 
T along the Euclidean time direction. Conservation of energy and angular momentum 
lead to two first order equations 



'''^ ihr , (5.5) 



^/gvPTJJWVh^ 

where the two conserved quantities are associated with the constants Uq and /, and the 
subscript indicates that the corresponding function is computed for u = Uq^ in which 
u{x) develops a minimum. We have used the notation /q = /(mq) and = h{uQ), and 
the prime denotes derivatives with respect to x. Solving these equations for x and 4> in 
terms of u we find 

^ = ^^^«"^i>'*/?S' (5.6, 

and 



F{u) ^ f{u) ( 1 - Ml ) _ (1 _ /^)/^ . (5.8) 



where we have defined 



h{u) 

Note that this function vanishes at m = uq, i.e. F{uo) = 0. If we place the quark at 

X = L/2, and the antiquark at a; = —L/2, the length of the Wilson loop is 



and 
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The total energy of the Wilson loop is divergent because it includes the infinite contri- 
bution from the W-bosons. When we subtract this contribution, the regularized energy 
of the quark-antiquark pair reads 



E,, ^ - Tdu J - - 1 r du^) , (5.11) 

where Uj^m is the minimum value of u allowed by the geometry. In specific examples we 
are supposed to solve for the auxiliary parameters Uq and I in terms of the separation 
distance L and the separation angle in the internal space A0. In practice this cannot 
be done explicitly for all values of the energy and length, but in some limited regions 
instead. In general, unless Z = we have that Acj) 7^ 0. If Z = then = and the 
angle remains constant. Then the function h{u) becomes irrelevant since it does not 
appear in the expressions for the length and the energy, and all our expressions reduce 
to those in [28]. In the specific examples that follow we will explicitly see that this is 
indeed the function encoding the deformation parameter 7. Therefore if / = our results 
will necessarily reduce to those for the undeformed theory. 

The heavy quark-antiquark potential as computed using the above formulas should 
obey the concavity condition 

dE - d'^E - 

§^>0, l§i<0. (5.12) 

stating that the force is always attractive and a non-increasing function of the separation 
distance. These conditions were proved quite generally in [41], and were investigated 
in detail for our examples in the case of zero deformation and angular parameters, i.e. 
for 7 = Z = 0, in [28]. For these values of the parameters and for certain trajectories 
the potential failed to obey (5.12) for all separation lengths. The failure is attributed 
to the fact that the trajectory approaches the singularity of the D3-brane background 
corresponding to the location of the branes and it is precisely this region that gives rise 
to the violation of the concavity condition beyond a certain length in the Wilson loop 
potential. Similar situations will also arise in this paper for non-zero deformation and 
angular parameters. 
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5.2 The conformal limit 



We will now analyze in detail some explicit configurations. In the conformal case all 
branes are located at the origin. We have 

,2 



g{u)^l, f{u)^u\ h{u)^^^^, (5.13) 
and Umin = 0. Also the function 

F{u) = {u^ - ul){e + ul{l - f)) . (5.14) 
Then we have for the length (we use 3.137(8 or 12) of [42])® 

L = R'vi^ui r 

[(2-f)E{k)-K{k)] , (5.15) 



where we have defined p = u^, and the modulus 



k = \lh^2- (5-16) 



For the angle we find (we use 3.131(8) and 3.137(8) of [42]) 

dp 



A(f) = V1+7V4 uo [ 



I ^ p{p - ui){p + ui{l - P)) 
= 2./l^W^^I==K{k) . (5.17) 

Therefore we see that the only effect of the deformation is a simple 7-dependent overall 
factor. The angle is a monotonously increasing function of the angular parameter I 
starting from zero. As I increases we reach its maximum value, 

A0„,ax = Vl + TT , for Z = 1 . (5.18) 

There arc special values of 7 for which the string has wound up n times the circle 

parameterized by the (j) angle. These are found by setting A^max — 27rn in (5.18). They 

are given by 7^ = 4(4n^ ~ 1)) n & Z. Of course the string can wound up for smaller 

''We will use the notation K.{k), E(fc) and n(n, k) for the complete integrals of the first, second and 
third kind, respectively. 



18 



values of / as well. It is not completely clear to us what the significance of these values 
for 7 is. For the energy we have (we use 3.141(12) and 3.141(18) of [42]) 
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1 

2^ 



dp 







7rV2^ 



{p-ul){p + ul{l-P)) VP. 
[K(fc) - (2 - /2)E(fc)] 



Uq 
TT 



(5.19) 



2i?2 [{2-f)-E{k)-}^{k)f l_ 

' L 



TT 



(2 - P)VT^ 

The Coulombic behaviour of the potential is characteristic of cases with conformal sym- 
metry since the only scale the enters in the various expressions is the quark- ant iquark 
separation distance. For / = this becomes the conformal result for J\f = 4 Yang-Mills 
[39, 40]. For larger angular parameter the result is still conformal but with an effective 
charge that is a monotonously decreasing function of I until it becomes zero for I — 1. 
This vanishing hmit corresponds to a BPS configuration [40]. 



5.3 The disc 

When vev's are turned one we expect that there should be some confining behaviour for 
the potential and from the gauge theory side this was studied in [22]. In the case of a 
brane distribution along a disc^ 

g(u) = 1 , f(u) = u\u' + rl) , h(u) = ""Y^^^f. 2 > (5-20) 



;i + 774)^2 + ro 



and Umi-n — 0. Also the function 



Fiu) = {u'-ul){u' + w), w^{ul + rl){l- ^y^tT/,^^Y\ ]>0. (5.21) 



It turns out that in computing the integrals we have to distinguish between two cases 
depending on which one of the two parameters w or Tq is larger. In fact we have that 



2 2 (l + 7V4)/2-l 

and similarly for the reversed inequality. Note that for small enough angular parameter 

we always have that w > r^, including the smallest value Uq — 0. It will be also convenient 

''Setting ro = (no brane distribution) and 7 = (no marginal deformation) reproduces the situation 
originally considered in [40]. The case Z = 7 = recovers the results in [28]. 
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to use the notation u'>(w<) to denote the larger (smaller) between the parameters w and 
Tq. Then we have for the length (we use 3.137(8) of [42]) 

dp 



where we have now defined the modulus 



n 



,k\-K{k) 



k 



w> — w< 



For the angle we find (we use 3.131(8) and 3.137(8) of [42]) 



A0 



Ih 



1/2 



/ 



1 + 7V4 + rl/p 



dp- 

§ \/ {p - ul){p + rl){p + w) 



2Ui 



1/2 



(l + f/4-roV«;>)K(fc) + 



n 



Also for the energy we have (we use 3.141(12) and 3.141(18) of [42]) 



,k 



Uq 
TT 



TT 



[(K^ + r2)K(^)-(K2 + ^^)E(A;)] . 



5.3.1 Generic behaviour 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



For Uq ^ To or, equivalently, at small distances, the vev can be ignored and the behaviour 
of the Wilson loop potential turns to be that of the conformal case (5.19).^ Towards 
larger distances the behaviour depends on the relation between the various parameters. 
If (1+7^/4)/^ < 1 then the potential becomes monotonically zero at a finite distance. This 
is the same as the behaviour of the Wilson loop in the undcformcd case corresponding 
to 7 = and for zero angular parameter / = found in [28]. If on the other hand 
(1 + 7^/4)^^ > 1 the potential, after it turns into positive values, reaches a maximum and 
then again approaches zero at a finite distance. In particular, using the properties of the 
^We should use the identity n(fc2, k) = E(fc)/(1 - k'^). 



20 



elliptic functions, we have that as ^ 



To V ^0 

A<^ ^ 7r/(l + ^^) , (5.27) 

Eg, - -[1-(1+7V4)/^]^, 
Hence we see that the energy goes to zero at a finite value of the length 

i^fin = ^VT^ , (5.28) 

resulting in a complete screening of charges at a finite distance. If we now use the 
expression for the length of the loop in (5.27) to solve for uq in terms of L, we find the 
vanishing behaviour 

Note that the way the zero energy is approached depends on the sign of 1 — (1 + 7^/4)/^, 
as was mentioned above. For a positive sign the above length Lgn is indeed the maximum 
distance between the quark and antiquark, beyond which there is no geodesic connecting 
them. In this regime the solution is described by two parallel strings with no interaction 
potential. The complete screening behaviour is qualitatively the same as that in [28], 
where the case with ^ — I — Q was analyzed in detail. In fact since the complete 
screening behaviour happens for — > 0, which is when the trajectory approaches the 
brane distribution, the continuous approximation breaks down and we should instead 
use the discrete version of the brane distribution. Then the potential becomes gradually 
Coulombic but with in (5.19) replaced by a much smaller value corresponding to the 
number of D3-branes located close to a particular center, i.e. i?^ ^T^gYM^ center i with 
-/^center <^ N . In that seuse the screening phenomenon still persists, but it is just made 
smoother. This conclusion assumes that N center ^ 1, so that the supergravity description 
is still valid. 

For a negative sign for 1 — (1 + 7^/4)/^ and depending on the strength of the parameter 
7 as compared with /, the energy is a single, or a triple-valued function of the distance 
(after a critical distance that can be computed numerically). We have depicted the 
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various behaviours in Figure 1. Beyond the maximum for the energy the force between 
the heavy quark and antiquark becomes repulsive and the concavity condition for the 
potential (5.12) is not obeyed. This region is not physical and we think that this behaviour 
is due to the large value of the deformations parameter 7 in relation to the trajectory 
approaching the singularity. Then of course the continuous approximation breaks down 
and the potential should reach a Coulombic behaviour as explained above. 



Energy 








Length 



Figure 1: The energy (5.26) as a function of the quark-antiquark distance L. As the value 
of 7 evolves, the shape of Eqq is modified and we go from the lower to the upper curve. 
For the upper two curves the result is trusted until the maximum energy is reached. In 
those cases the maximum length Lmax < -^fin- In the lower curve Lmax = L^n- 



5.3.2 Confining behaviour 

In order to avoid the apparently unphysical region after the peak of the potential we 
should cut off the region beyond this critical point. A careful analysis shows that this 
amounts to avoiding the region in the deep IR by essentially restricting the energy pa- 
rameter to Uq ^ ^o/t- In the limit of large 7, the potential and length are given by 
(5.23) and (5.26) but with w — {1 — l'^){u1 + Tq) which arises from the definition (5.21) 
in the limit 7 — > 00. We easily find that the potential reaches a constant positive value 
at a finite value of the separation distance, thus again resulting in a complete screening 
of charges. They are given by 

£;.,,fin = - [K(0 - E(0] , Lfin = ^ . (5.30) 

TT ro 

The positive value of the maximum energy depends on the angular parameter / in such 
a way that it grows when I approaches 1 as Eqq^^^ ~ — |^ln(l — Z^) — > 00, so that in 
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practice the potential remains confining. This aUows to consider the hmit of a very large 
deformation parameter and simultaneously take the angular parameter to one, that is 



7 > 1 and 1^1 . (5.31) 
This introduces a hierarchy of widely separated scales and in particular note that 

Next we will restrict to the potential corresponding to the corner of the parameter space 
(5.31) and further impose the condition 

? « no « -yfi= . (5.33) 
7 V 1 — t 

Since the parameter uq plays essentially the role of the probe energy, this means that we 
will not probe with extremely high energies in the UV, i.e. extremely higher than the 
vev value or, equivalently, the separation length will not become extremely small. This 
cuts off the conformal region. Similarly, by never probing with energies in the extreme 
IR we decouple the region after the maximum of the potential. We find that in the limit 
(5.31) the length (5.23) becomes 

- 2 E{k) - k'^K{k) 



ro kk' 
where now the modulus simplifies to 



(5.34) 



k^ J ^ ° ^ , with k' = Vl-A;2 (5.35) 
V «o + ^0 

and we have defined the finite length 

L = ^ . (5.36) 

For the angle (5.25) we find 

A0 = 7 K(A;) , (5.37) 
and for the energy (5.26) we obtain 

E„^lim^, (5.38) 

The potential obeys now the concavity conditions (5.12), and its behaviour is depicted 
in Figure 2. 
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Figure 2: The energy (5.38) as a function of the quark-antiquark separation. For large 
values of the distance we find the logarithmic confining dependence (5.43). For smaller 
separations the behaviour is linear, (5.40). 

Consider the large values 

rn - 1 - a/1 - P 

uq > ro (and uq < ) ^ L -C — (and L > ) , (5.39) 

yl — P fo To 

where the inequalities inside the parentheses emphasize that energies cannot become 
extremely large and correspondingly the lengths cannot be too short. Otherwise the 
limit (5.31) that resulted into (5.34)-(5.38) is not self-consistent in a mathematical and 
a physical sense. Then we obtain for the energy a linear behaviour as a function of the 
separation distance, 

E,,c^l^^^L^O. (5.40) 
4uo 27r 

For the angle we obtain the large limiting value 

A0 = 7 I . (5.41) 

Let us now consider the opposite limit of low values 

rn - 1 - 7 

uo -C ro (and Mo > — ) <^ L > — (and L < — ) , (5.42) 
7 ^'o ro 

where, similarly to the previous case, the inequalities in the parentheses imply that the 
region in the deep IR is avoided for the hmit (5.31) to be self-consistent. For the energy 
we find 

E,. ~ In f ^ V !l In (rol) ^ oo . (5.43) 



24 



Hence we have a confining behaviour for the heavy quark- ant iquark potential, ahhough 
the dependence on the separation distance is not hnear but logarithmic. For the angle 
we obtain again a logarithmic behaviour 

~ 7 In(roL) . (5.44) 

It is worth commenting that a logarithmic form for a confining potential instead of a linear 
one is not strange to high energy physics. Long ago it was shown that a logarithmic 
potential is unique in producing spectra with energy level spacing independent of the 
particle mass [43] , an idea that arose from related observations at particle spectra at the 
time. In the same work the logarithmic potential was used in relation to quarkornium 
spectra for various families of heavy qq bound states (for related reviews on the use 
of Schrodinger equation in heavy quark systems see [44]). This phenomenological way 
of introducing a logarithmic potential in Particle Physics had no theoretical support 
within perturbative approaches to QCD. More recently, in the context of potential non- 
relativistic QCD, logarithmic terms arose, in calculations valid for length scales less than 
1/AqcD) as modifications of the effective charge in the Coulombic potential (see [45] and 
references therein). Given our results, it will be interesting to investigate the issue of a 
logarithmic confining behaviour from a gauge theoretical view point. 



5.4 The sphere 



This case follows from the expressions above if we let — > —r^ and take into account 
that now limin = '"o- Therefore we have 



f{u)=u\u'-4) , h{u) 



(1+ 774)1^2 -rg 



(1+774)/^^^^ 



(5.45) 



(5.46) 



g{u) = 1 , 
Also the function 

F{u) = iu^-ul){u^ + w), (l + 774)«g 

It turns out that in this case we always have, although w is not a strictly positive constant, 
that w > — rg, and therefore we do not have to distinguish between two different cases 
depending on the values of w and Tq, as in the disc case. 

Then we have for the length (we use 3.137(8) of [42]) 

dp 



J u; 



w 



ul + w 



n 



w 



,k)-K{k) 



(5.47) 
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where we have defined p = u^, and the modulus 



w + rQ 
w + ul ' 



(5.48) 



For the angle we find (we use 3.131(8) of [42]) 



1/2 



dp 



l+7V4-r2/p 



V'(p-^o)(p-^o)(P + ^) 



2lk 



1/2 



w 



w 



w 



,k 



(5.49) 



For the regularized energy we now have (we use 3.141(12) of [42]) 



QQ 



1 



{p-ul){p + w) ^ 



Up - rp 

TT 



^== [{ul - rl)V.{k) - {ul + wMk)] + 5 



(5.50) 



As before, for short distances we may ignore the vev and the behaviour of the Wilson loop 
remains that of the conformal case. For larger values of tq we may easily show, using the 
properties of the elliptic functions, that the length goes again to zero and therefore the 
energy is a double valued function of the separation length. As before, if we consider the 
limit of large 7 together with I ^ 1 we will cut off the conformal part of the potential. 
In these limiting case we obtain for the length and the energy the simplified expressions 



and 
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^=^|[K(fc)-E(fc)] 



^(^K(.)-1e(.)^1 



(5.51) 



(5.52) 



where now the modulus is k = tq/uq and L is defined as in (5.36). Using properties of 
the elliptic functions one easily verifies that the separation L goes to zero both for /c ^ 
and A; — > 1 and therefore there is no confining behaviour. We note that there is such 
behaviour (linear) for the — ^ trajectory which however, as remarked in footnote 5, is 
not sensitive to the deformation. The energy is a doubled valued function of the length. 
The angular distance now becomes 



A0 = ^k'K{k) . 



(5.53) 
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6 The wave equation 



In this section we study the massless scalar field equation 

= ^^e-'^^M^^e-^^G^^^AT* = , (6.1) 
y—G 

in the deformed SO{A) x SO{2) background. Quite generally, for the background metric 
(2.24) and the dilaton (2.26), it can be shown that the above equation can be written as 
a deformation of the corresponding equation for the undeformed background, 

= 0^=0* + 7^ Ay didj'^ = , (6.2) 

where we have defined 

_ Z1Z2 + 2^3 + ^2^3 r Z1Z2Z3 

Xij = dij - — — . (6.3) 

Let us now study each of the pieces entering this differential equation. 



6.1 The undeformed case 

First we will examine the Laplacian for the undeformed case with 7 = 0. This case is 
on its own already very interesting because it describes the spectrum of dilaton, trans- 
verse graviton and gauge field fiuctuations in the undeformed background.^ Within the 
AdS/CFT correspondence these fluctuations correspond to gauge theory operators. In 
that respect the s-wave case has been explicitly first solved in [27, 28]. In that case 
the description can be made also using the five-dimensional gauged supergravity arising 
from reduction of type-IIB supergravity on S^. Here, we will see that we can go beyond 
gauged supergravity and explicitly solve this equation in all generality. We will shall see 
below that, turning on the deformation will not affect dramatically most of our analysis 
of this subsection. 

For the background corresponding to the disc distribution we make the ansatz 

* = (^^*''^^^^' ^' '^i' ^' '^2' '/•s) • (6.4) 

^These spectra are in fact degenerate and the corresponding fluctuations are interrelated and can 
be written in terms of the dilaton fluctuations. These have been exposed in a series of papers, for the 
dilaton-gravity in [46] and for the dilaton-gauge fields in [47, 48] . This interrelation is essentially due to 
the fact that all fluctuations correspond to fields that belong to the same Af = 4 supergravity multiplet 
[47] . It would be very interesting to see whether such a coincidence of spectra persists in the marginally 
deformed Af = 1 case as well. 
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Obtaining the corresponding differential equations below for the sphere or the conformal 
cases simply amounts to analytically continuing the parameter Tq to —Tq or setting it to 
zero, respectively. Then after some algebra we find the following second order differential 
equation for the Laplace operator acting on the amplitude 



□ 0^ ~ 



+ 



1 rl 



R^W 1^3,2 2^ o 1 

H smt^cos'^ty 



d\ + f ^ + % Ass 



^sin^^ r2 + r2y '^^ 
where the mass eigenvalue is defined as 



sin 9 cos^ 9de 



(6.5) 



(6.6) 



and where we have omitted an overall coefficient and the factor from the exponential. 
Obviously the Laplace equation corresponding to the operator in (6.5) admits solutions 
via the separation of variables method. Let 

1 



^±(r,^,0i,V',02,03) 



neZ . 



For the ^53 piece we get the standard eigenvalue equation on 

A53*53 = -/(/ + 2)^53 , / = 0,1,... , 

or, explicitly. 



sin 21/; dijj \ ^ dip J 
The normalized solution is given by 



sin 2^::^ ) + (/(/ + 2) + — f-^ + ) = . 



cos^ sin 



(6.7) 



(6.8) 



(6.9) 



^S3,;,n„n3(^,02,</>3) = A;,„,,„,e'("^^^+"3<^3) (^i^ ^) l"3l (^os ^) l"^lpf (cOS 2^) , 



I — 1^2] — Irial — 2k , k — 0,1, ... , 
where the normalization constant is 



l,ni,m 



/ + ir(i/ + i 


n2 


+ i 






n2 


1 

2 




+ 1) 


TT r(i/ + i 




1 

2 


n2 


+ i)r(i/-| 




+ \ 


n2 


+ 1) 



(6.10) 



(6.11) 



The mass eigenvalue has degeneracy 
d3,i ■■ 



— ^2k+\n2\ + \n3\,l - S2k+p,l ^ ^|n2| + |n3|,p 

fe=0 n2,n3=— oo k,p=0 n2,n3=— oo 

oo 

= ^ {Ap + Spfi)S2k+p,i 

k,p=0 
= (^ + 1)% 



(6.12) 
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which is the correct expression. Using the above we obtain, for the other factors, 



sin 6 cos^ 6 d6 

and 



,in9cos=0!?^+fB-4^-^Vm = (6.14) 



iM^R" - 1(1 + 2)rl)r + -^^n" - Er^ 

+ Vq 



= , (6.15) 



where E is the separation of variables constant. The measure for the resulting Hilbert 
space is 

v^G" = d'^x{drr) {d9 sin 9 cos^ 9) {d^j sin ^ cos ^) . (6.16) 

It is quite remarkable that the equations for the angular part, and in particular (6.14), 
do not depend on the parameter tq characterizing the vev of the scalar fields. Hence 
we expect that E should be quantized and written as E — j{j + A), j — 0,1, .. . as 
appropriate for the Laplacian on the undeformed S^. This is not immediately apparent, 
but it works precisely like that as we shall prove. Writing 

V'^'^ = (1-x)I"I/'(1+x)'/2F(x) , x = cos2e, \x\<l, (6.17) 

we obtain an equation for F{x) which is the Jacobi differential equation with a = |n| 
and P — I + 1 in the standard notation. It has the Jacobi polynomials as a complete 
set of orthonormal solutions provided that the parameter E is quantized as 

Em,i,n^ {l+\n\+2m){l + \n\+2m + A) , m = 0, 1, . . . . (6.18) 

general eigenstate of the Laplace operator on the unit n-sphere has energy eigenvalue and degen- 
eracy given by (see, for instance, [49]) 

End = j{j + n-l) , j = 0, 1, . . . , 

(2i + n-l)(i + n-2)! 

(n-l)!j! • ^^^^> 

^^More precisely, if we change variable as z = (1 + cos20)/2, then the differential equation for ^^^^ 
transforms into a second order ordinary differential equation of the Fuchsian type with three regular 
singularities at z = 0, 1 and oo. Such an equation can be transformed into the canonical form of a 
hypergeometric equation. The latter is transformed to the Jacobi equation by a simple variable change 
X = 2z — 1{— cos29 as in (6.17)) which has polynomial solutions provided the parameters a, h, c in the 
standard notation are related to, at least, an integer. The same reasoning applies for the differential 
equation for ^'^^ that is considered separately for the conformal, the disc and the sphere cases below. 
We explain in section 6.3 that in the deformed case the nature of the singularity at 2: = oo changes from 
a regular to an irregular one and as a result we cannot solve the corresponding differential equation by 
elementary methods. 
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Indeed this takes the form j{j + 4), with j = I + \n\ + 2m, as advertised. Moreover, 
the degeneracy of this state works out as it should be. Indeed, note that in general the 
degeneracy is given by 



oo oo 



m,l=0 n=— oo m,p=0 



oo oo 



1=0 n=— oo 

= ^ ^ (p + 1) (2p2 + 4p + 3)S,+2m,j (6.19) 

m,p=0 

- ^0'+i)a+2Hi+3), 

which is the correct expression (see (6.13) with n — 5). Then 

<ln(^) = 5.n,/,nsinl"l ^cos' ep(l"l''+^)(cos 2^) , (6.20) 

where the normalization constant obeys 

^9 / , , ,,r(m+ l)r(m + Z + Inl + 2) 

52 . = 2 2m + 2+ + n , \ ^ — -( . 6.21 

^ ' '^r(m+ |n| + l)r(m + Z + 2) ^ ' 

The integration measure that determined the normalization constant is d^sin^cos^^, 
with < ^ < 7r/2. 

Let us now present the solution to the remaining radial equation. We will have to 
distinguish between the conformal, the disc and the sphere cases. 

6.1.1 The conformal limit 

In this case we have to solve equation (6.15) after setting the parameter ro = 0. This 
can be written as a Schrodinger equation by letting 

^(2) ^ ^3/2^/ N ^ z=- , z>0 . (6.22) 

r 

Then ^ satisfies the standard Schrodinger equation with potential 

and eigenvalue M'^R'^. This potential is positive definite and vanishes for large values of 
z. Hence we expect a continuous spectrum with no mass gap. The explicit solution is 

V'(') = v^^V,+2(Mz) , J = 0, 1, ... , (6.24) 

where j is the integer that parametrizes the eigenvalue in (6.18) as we have mentioned 
above. 

30 



6.1.2 The disc 



In the previous section we have found a vanishing behaviour beyond a maximal length 
for the interaction potential between a heavy quark and an antiquark. A maximal length 
implies the existence of a set of small energy scales that we can not reach. Therefore 
the gauge theory must generate a mass gap at strong coupling. We will now prove the 
existence of this mass gap from the solution to equation (6.15). We note first that (6.15) 
can be written again as a Schrodinger equation. Indeed let 

^(2)^sinh^^ sinh;2=^, z>Q. (6.25) 

Then ^ satisfies the standard Schrodinger equation with potential 

V{z) = (/ + 1) -2— + ■ , 2 ' (6-26) 

cosh z smh z 

and eigenvalue M^R^ lr\. This belongs to the family of Poschl-Teller potentials in quan- 
tum mechanics of type II. The potential decreases monotonically from arbitrarily large 
positive values, where it behaves as in the conformal case in (6.23), to the constant (/ + 1)^ 
as z ranges from to oo. Hence we expect a continuous spectrum with a mass gap given 

by 

^gap,. = (^ + l)§, (6-27) 

which is degenerate according to (6.12). In order to explicitly solve equation (6.15) we 
perform the change of variables 

^(2) ^ ^-1/2+^2(1 _ Qm^\n\l2^2p^^^ ^ ^ ^ ^_ Q < ( < 1 , (6.28) 

+ Tq cosh z 

where 5 is a constant to be determined. After some algebra we obtain a standard hyper- 
geometric equation for F(C), whose general solution leads to 

V'S,„(r) = AAC-^/2(l-C)"^+l"l/'+Me*T'/Vi(a,a,l+i<^,C) + c.c.) , 

= AAC-'/'(l-C)""'"'"'^'(e'T^/'F(^',^',l+^ + ^5,C) + c.c.) (6.29) 

where M is an overall normalization constant and 

a= l(|n| +2m + 3 + i5) , 6 = 1 + i5 - a , b = Jw^R^/rl - {l + lf . (6.30) 

The phase </? is computed by demanding that the solution is regular at C = 1 (equivalently 
at r — > oo). We find that 
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As in [27, 28], where the s-wave case with I = n = was studied in detail, normahzabiUty 
of the solution in the Dirac sense requires that the parameter 6 is real and therefore the 
spectrum is continuous, with the mass gap (6.27). 



6.1.3 The sphere 

As in the previous cases eq. (6.15) (with Tq —>■ —Tq) can also be written as a Schrodinger 
equation. Now let 

• 3/2 

(2) ^ sm ^ sinz^- , < ^ < 7r/2 . (6.32) 

cos^/^ z r 

Then ^ satisfies the standard Schrodinger equation with potential 

V{z)^-{1 + 1) + — + — '- , 6.33 

cos^ z sm z 

and eigenvalue M'^R^/tq. This belongs to the family of Poschl-Teller potentials in quan- 
tum mechanics of type I. To find the explicit solution we perform a change variables 
again through 

^(2) = (1 _ ^)m+;/2+|n|/2(-^ ^ C)'"'^'^(C) , C = 1 " 2^ , |C| < 1 ■ (6-34) 



Then we obtain the standard Jacobi differential equation, with a = 2m + I + n — 2 and 
P — n. The corresponding normalized solutions are 

c,,.w - c,^,™,. (?)'"'""'"'" (i-^)'"'%r*'*'^^^'""(i-2.s/.^) . 

A; = 0,1,... , (6.35) 
where the normalization constant obeys 

- 2(2h I 2m III 2\n\ I 3) + mk + 2m + I + 2\n\ + 3) 1 

C,,m,n - 2[2k + 2m + / + 2|n| + 3) ^ 2m + / + H + 3)r(A; + H + 1) ^ ' ^^"^^^ 

The mass eigenvalue is 



Km,i,n = 4(/c + m + |n| + l)(/c + m + |n| + / + 2)-^ . (6.37) 

The measure that determined the normalization constant is rfrr, with ro < r < oo. This 
mass eigenvalue is of the form = 4(j + l)(j + / + 2)rl/R'^, with j = k -\- m -\-\n\. Its 
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degeneracy follows from a computation similar to that in (6.19), 

h+\n\,p 



oo oo 



<^sphere,ij — ^k+m+\n\,i — (^ + 1)^ ^m+pj 

m,fc=On=— oo m,p=0 



oo oo 



.A;=0 n=— oo 



= (Z + 1)^ ^ (2p + l)5^+,,,. (6.38) 

m,p=0 

= (/+l)^(j + l)^ 



6.2 The deformed case 

In this case we have to take into account the second term in the right hand side of (6.2). 
In general this term destroys the 5'0(4) spherical symmetry, unless we consider solutions 
that are independent of the angles 02,3 of 5"^, as parametrized in (3.9). We will then 
consider the Laplace operator within the ansatz 

*=7^e*'*^(^'^''^i'^) ' (6-39) 

which is a restricted version of (6.4). We obtain now an expression similar to (6.5), with 
an additional term inside the bracket. In particular. 



+ 7^cos^^a^, 



*± , (6.40) 



R^yjr"^ -\- rl cos^ 9 

where the eUipsis denotes the terms inside the bracket in (6.5). The measure in the Hilbert 
space is given generally by V— Ge~^*G**. This is given, based on general grounds, by 
the right hand side of (6.16), but it can be also verified by an explicit calculation. 

Proceeding as in the undeformed case we make for the amplitude the ansatz^^ 

*x(r,^, 01,-0) = -^e^"'^i*53(^/')^/^(^)(^)^/^(2)(r) , n e Z . (6.41) 
V 27r 



For ^53 we obtain the same equation as in (6.8) but with no 02,3 for the solution. Hence, 
from (6.10) we have that Z = 0, 2, . . . so that we make the replacement I — > 21. Now 

^^Since e~'^^\fG is, in general, invariant under T-duality and G** is not affected by the T-dualities we 

performed. 

-'^■^If we make a general ansatz of the form ^ e™''^'-0'^^^ (-0, 9)ip'-^^r), where rii G Z, the modification 
of the differential equation depends on the matrix Ajj in (6.2) via a term of the form —^^XijUiUj. For 
general integers rij this ansatz will break the spherical 50(4) symmetry of our background and there 
is no longer factorization of the V and 6 dependencies. However, if 722 = the spherical symmetry is 
preservcid. Then using (3.11) we may easily show that the expressions below are still valid provided we 
replace by (n — 71.2)^. 
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the eigenvalues are of the form 4/(/ + 1) and the normaUzed solution is given in terms of 
Legendre polynomials as 

/ IX 1/2 

^SM = 2il + -\ P;(cos2^) , / = 0, 1, . . . . (6.42) 

For the radial factor V'^^^ we obtain the same equations as in (6.15) with the replacement 
/ 21. The only major modification is in the equation for the angular part which 
now gets an additional term and becomes 



1 



d 



sin 6 cos^ 6 dO 



sin d cos d 



dAlP^ 

de 



E-n' 



1 ^2 2 A 4Z(i + l) 
^^+7^cos^g - ' 
sm 9 I cos^ 9 



This equation can also be written as a Schrodinger equation. Indeed let 

sinV2^ cos3/2^ ■ 
Then ^ satisfies the standard Schrodinger equation with potential 



V^(^)= 0.(6.43) 



(6.44) 



V{9) = -4 + 



n 



1/4 4/(/ + 1) + 3/4 



+ 



+ n^7^ cos^ 9 



(6.45) 



sin^ 9 ' cos^ 9 

and eigenvalue E. When n = (or n = ^2 = 713, see footnote 13), the results obtained for 
the undeformed case carry over unchanged. For n 7^ we could not solve this equation by 
elementary means for reasons related to the change in nature of the singularity at infinity, 
as described in footnote 11. Nevertheless, we may resort to perturbation theory which 
is valid for small values of the effective parameter n^7^, and to an asymptotic expansion 
for large values of the same parameter. We will return to the generic case later. 



6.2.1 The case with 7^71^ < 1 

In the hmit 7^n^ <S 1 we can treat the last term in the potential (6.45) as small. The 
corresponding shift in the energy eigenvalue can then be found from conventional per- 
turbation theory. Using the wavefunction (6.20) we get 

5E^,i,n = ln^'Bl,^,j\x{l-x)\-\l + xy^'[Pj:\^^^\x)r >0. (6.46) 
For I — the integral can be computed. We found that 



SE n — 2^+1"! 



(m 


+ 


l)(m + 


n 


+ 1) 


(2m + 


n 


+ 3)(2m + 


n 


+ 1) 



(6.47) 



The corresponding change in the eigenfunctions ^^^^ is also computable using perturba- 
tion theory, but it will not be presented here. The change in the values of E^^i^n a-re 
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affecting the differential equation for the radial factor 'tp^'^^ according to (6.15). It can 
be easily seen, for instance from the form of the potentials for the effective Schrodinger 
problems (6.23), (6.26) and (6.33), that the effect due to the shift 5-Em,«.n can be absorbed 
by an effective shift in the quantum number m. This can be computed by varying (6.18) 
with respect to m and identifying the result with the perturbative result (6.46). In this 
way we find 

Sm = . (6.48) 

4(/ + |n| +2m + 2) ^ ^ 

For the conformal and the disc cases the continuum mass spectrum for remains 
unchanged, except of course for the density of states which depends on the wavefunction. 
For the sphere case the spectrum changes as 

2fc + 2m + 2|n| + / + 3 „ 

6.2.2 The case with j'^n'^ > 1 

In this case for generic values of the angle 6, i.e. not near 6 = 0,7r/2, the potential goes 
to +00 and the wave function has the form of plane waves of very high energy. However, 
there are special states that are actually quasi-localized or even localized near ^ = and 
9 — 7^/2, respectively. 

• Region near ^ = 0: Let the change of variables 

. (6.50) 

77-7 

Then, within our limit the variable z becomes non-compact with 2; > and the problem 
becomes equivalent to a Schrodinger problem with potential 

V{z) = VLzlll + 1 , ^ > , (6.51) 



so that the spectrum is continuous with mass gap. The solution is 

E 

with the measure induced by (6.16) being dzz. Note also that in our hmit i)^^^ ~ z"^!"^^. 
This solution is quasi-locahzed near 2; = (for ^ ~ ^) in the sense that it dies off 
oscillating with a power law behaviour away from it. 

• Region near 6 — 'k/2: Let the change of variables 

(6.53) 
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Then, as before, within our limit the variable z becomes non-compact with 2; > and 
the problem becomes equivalent to a Schrodinger problem with potential 

^(,).!i(i±4±M + .2_ .>o (6.54) 

and rescaled energy parameter E — The essential difference with the previous case is 
that the potential gets a term corresponding to a harmonic oscillator. The Schrodinger 
equation can be transformed into a confluent hypergeometric equation. The solution 
regular for 2; = is then given by 

^(1) ^ z^^e-''^^F{l + 1 - E/A, 21 + 2, z^) . (6.55) 

The measure induced by (6.16) is dzz^. Note also that in our limit ~ z"^/^^. For 
generic values of E, the confluent hypergeometric function appearing in the solution 

2 

diverges exponentially as , as 2; — > 00. In order to avoid this behaviour at infinity 
we demand polynomial solutions which implies a quantization for the eigenvalue E, as 
expected from the shape of the potential. Writing the result for the original energy 
constant E we have that 

~ 4(m + Z + 1) ^7 , m = 0, 1,... as ^7 > 1 , (6.56) 

When 7^1 this formula is vahd for n > 1. Obviously, the solution is locahzed near 
z — {) (for I — ^ ~ in the sense that it dies off exponentially away from it. 

6.3 Relation to the Heun equation and the Inozemtsev model 

We turn now to a general study of the equation (6.43) for the angular function ■0(^)(^) 
in the deformed case. The substitution 

V'^'^ = (1 - ^)I"I/'^'F(^) , z = cos' e , (6.57) 

gives the following equation for the function F{z) 
d'F ^ r 2(1 + 1) ^ \n\ + l \dF ^ n^'z-E+{2l + \n\){2l + \n\+A) ^ ^ ^^^^^^ 



dz^ \ z z — 1 J dz 4:z(z — 1) 

This is a second order ordinary differential equation with two regular singular points 
at 2; = and at 2; = 1, together with an irregular singularity ai z = 00. Hence it is 
expected to correspond to a conffucnt form of the Heun differential equation (see, for 
instance, [50]). The latter is the standard form of a Fuchsian differential equation with 
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four regular singularities. Similar to the case of the confluent hypergeometric equation, 
under a confluence process two of the singularities are made to coincide resulting into an 
irregular one. 

There is an appealing relation of the Heun differential equation to the integrable BCi 
Inozemtsev model. The BCi Inozemtsev system is a one-particle quantum mechanical 
model in one dimension with potential 



i=0 

where p{z) is the Weierstrass elliptic function, which is doubly periodic in the z variable 
with half-periods a; and a;'. The Weierstrass elliptic function obeys 



where the constants e^'s obey p(a;) = ei, p{u! + 00') — 62 and p(a;') = 63. In case they 
are in a straight line in the complex plane, 62 is assumed to lie between ei and 63. In 
addition, they sum up to 0. In the notation of (6.59) ujq — 0, cui — cu, L02 — uj' and 
cus — cu + ijj' . This model was shown to be completely integrable in [51] and it belongs 
to the class of -BCat integrable quantum N-particle systems with symmetry. It is a 
generalization of the -BCat Calogero-Moser-Sutherland model. 

The expression of the Heun equation in terms of elliptic functions was essentially 
known by Darboux. The explicit relation of the BC\ Inozemtsev model to the Heun 
equation can be found, for instance, in [50, 52, 53]. The potential (6.59) was also studied 
in [54] , in the context of the theory of eUiptic solitons. It was shown that if the above 
coupling constants Zj, with i = 0, 1, 2, 3, are all integers, then the potential (6.59) has the 
finite gap property. 

We will now consider the so-called trigonometric limit of the Schrodinger equation 
corresponding to (6.59). In the standard version of this limit we end up with the trigono- 
metric Poschl-Teller potential, which in our case is the potential (6.45), but with 7 = 0. 
Therefore, if such a limit is taken, we do not capture the deformation. Here we will 
take instead a type of trigonometric limit which can take into account the effect of the 
deformation. In the mathematics literature this was considered before in various forms 
(see [55, 56]). In general, in the trigonometric limit, we take the limit —> 63. Then 
one of the half-periods becomes imaginary infinity, say uj' too. In order to investigate 
this limit it is better to consider the series representation of the Weierstrass function in 
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(6.59) 



{p'f = 4(p - ei)(p - e2)(p - 63) , 



(6.60) 
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terms of powers of the parameter q = e*'^^, where r is the modular parameter defined as 
r = uj^/ui. We take the half-periods to be cji = |, c<j2 = and c<j3 = —Ui — uji = |. 

Then using standard expansions for the Weierstrass function and keeping only the rele- 
vant to our discussion terms, we find 



2 2 

p{x) = + -4— + 167r' sin^ tix + O(g^) , 
6 sm TTX 



2 2 

p(a; + a;i) = + + IGtt^ cos^ tto; + O(g^) , 

3 COS"' TTX 



p{x + U2) = h87r^cos27rx g + 87r^(l - 2cos47rx) r + , (6.61) 

3 

p(x + U3) = — - - 87r^ COS 27rx 5 + 87r^(l - 2 cos47rx) q'^ + 0{q^) . 



Next we scale and shift the energy as 

E 

and let 



^'(^E-^J2j^{k + l)j (6.62) 



where c, d are real constants. Then the limit g ^ exists even though the parameters I2 
and I3 in the potential both go to infinity. If we define the angular variable 9 = —n/2+T:x, 
taking values in the range 9 e [— 7r/2,7r/2] we obtain a one-dimensional Schrodinger 
equation with potential 

V{9) = + + M!l±i) + 2c\l - 2cos4^) - 16(icos2^ . (6.64) 
cos^6' sin^6' 

We call this the generalized trigonometric limit, in distinction to the ordinary one for 
which c = d = 0. The potential (6.45) corresponds, up to an additive constant, to the 
particular case c = 0, d = —^n^l"^, = 2/ + | and h = n — \. Also the range of 9 
is restricted to ^ e [0, 7r/2] (this can be done by requiring vanishing eigenfunctions for 
^ = 0). 

The advantage of having a precise relation of (6.45) with the Heun equation and 
the BC\ one-particle quantum Inozemtsev integrable model is the fact that techniques 

have been developed based on the Bethe ansatz in order to determine the eigenvalue 
problem for the general potential (6.59). In particular, in [57] an explicit solution was 
presented. However, it still involves satisfying / -|- 4, with / = Y^i=o ^ii complicated 
relations involving Theta-functions for I auxiliary parameters. Solving this problem in 
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general seems impossible. In [57] it was shown how to solve this problem in the ordinary 
trigonometric limit employing the Bethe ansatz. The result reduces to Jacobi polynomials 
as in our case, in a presumably equivalent solution. We expect that progress can be made 
using the Bethe ansatz method even in the generalized trigonometric limit. 

7 Conclusions 

In this paper we have constructed type-IIB supergravity duals to the Coulomb branch of a 
class of exactly marginal deformations of A/" = 4 supersymmetric Yang-Mills theory. The 
solutions can be derived by applying a sequence of T-dualities and coordinate shifts to 
multicenter supergravity backgrounds, similarly to the conformal case [6] . On the gravity 
side of the AdS/CFT correspondence the marginal perturbation leads to a deformation of 
the space transversal to the worldvolume of the branes. We have probed the geometry of 
the marginal deformation by computing the expectation value of the Wilson loop operator 
from a string extending along a great circle of this deformed transversal space. The cases 
that we have considered in detail correspond to a continuous distribution of D3-branes 
on a disc or on the surface of a three-dimensional sphere, and preserve an 5*0(4) x 
5*0(2) global symmetry group. In the conformal limit, where the Higgs vev remains 
small, we observed the usual Coulomb behaviour for the quark-antiquark potential. The 
background corresponding to a distribution of branes on a disc is quite rich, and contains 
several regimes according to the relation between the various parameters. In particular 
we have found situations where the quark-antiquark interaction is completely screened 
as in [28], or where a confining behaviour arises, with a linear or logarithmic potential 
according to the ratio of the quark-antiquark separation distance to the Higgs vev scale. 
It would be of great interest to explore this novel behaviour on the field theory side of the 
correspondence and in particular to explore the origin of this logarithmic dependence. 

We have also described the spectra of massless excitations by solving the Laplace 
equation in the deformed 5*0(4) x 5*0(2) background. In order to solve the correspond- 
ing differential equations we have transformed them into Schrodinger equations with a 
potential determined by the geometry of the supergravity background. We have per- 
formed a detailed analysis of the equation in the J\f — A limit, with no deformation, 
for the conformal, the disc and the sphere distributions. It is quite remarkable that it 
is possible to explicitly solve this problem beyond the s-wave approximation that has 
been considered before [27, 28]. It will be interesting to extend the analysis to the 
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full class of models corresponding to the gravity-scalar sector (with scalars in the coset 
150(6, ]R)/ S0{6)) of five-dimensional gauged supergravity that admit a classification via 
a correspondence with algebraic curves [29]. In the deformed background we have relied 
on perturbation theory to find a solution for small values of the deformation parameter, 
or on an asymptotic expansion in the opposite limit. Generic values of the deformation 
lead to a confluent form of the Heun differential equation, which is known to be related 
to the Inozemtsev integrable system. This relation provides a tool to find solutions to 
the differential equation through the Bethe ansatz method. In this article we have intro- 
duced a modification of the usual trigonometric limit to give account of the deformation 
parameter. A complete attempt to investigate the spectral problem in the generalized 
trigonometric limit that we have performed is definitely worth the effort. 
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